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1.  Introduction 


The  unique  properties  of  ferroelectric  materials  have  led  to  numerous  studies  aimed  at  incorpora¬ 
ting  ferroelectrics  into  memory  devices.  One  particular  application  is  the  creation  of  non-volatile 
ferroelectric  random  access  memory  (NVFRAM).  A  ferroelectric  material  can  maintain  an 
induced  polarization  without  an  externally  applied  voltage,  thus  removing  the  need  for  a  back-up 
disk  and  making  the  memory  non-volatile  (i).  The  ever-increasing  demand  for  smaller  memory 
devices  has  led  to  increased  research  in  the  field  of  ferroelectric  nanostructures.  Ferroelectric 
nanodots  offer  a  potential  increase  in  the  storage  density  of  NVFRAM  from  0.2  gigabit  per  square 
inch  (Gbit/in.  )  to  60  terabits  per  square  inch,  which  is  also  four  orders  of  magnitude  greater  than 
the  1 -Gbit/in.2  capacity  of  current  magnetic  recording  devices  (2).  The  improved  capacities  for 
memory  have  inspired  much  interest  in  the  properties  and  behavior  of  ferroelectric  nanostructures. 
One  area  of  interest  is  the  toroidal  ordering  of  dipoles  in  ferroelectric  nanodots.  This  is  a  relatively 
unexplored  field  that  appears  promising  yet  requires  extensive  research  before  commercialization 
can  be  seriously  considered. 


2.  Background  on  Ferroelectrics 


Ferroelectric  materials  are  in  a  unique  category  of  materials  defined  by  the  appearance  of  spon¬ 
taneous  polarization.  Ferroelectric  crystals  are  a  class  of  dielectric  materials,  which  means  they 
are  semiconductors  that  have  insulating  and  conductive  properties.  They  may  sustain  internal 
electric  dipole  moments  without  an  externally  applied  electric  field.  By  contrast,  dielectrics 
generally  require  an  electric  field  to  sustain  a  polarization.  Ferroelectric  properties  also  make 
them  attractive  for  numerous  applications  in  multifunctional  electronic  devices.  Such  properties 
include  high  dielectric  constants,  the  pyroelectric  effect,  the  electro-optical  effect,  and  the 
piezoelectric  effect. 

Ferroelectric  crystals  become  polarized  when  the  positively  charged  central  atom  shifts  in  the 
opposite  direction  of  the  negatively  charged  oxygen  atoms.  Figure  1  illustrates  how  a  non¬ 
polarized  ferroelectric  crystal  becomes  polarized.  This  phase  transition  only  occurs  below  a 
critical  temperature  referred  to  as  the  Curie  temperature  or  Curie  point.  The  value  of  the  Curie 
temperature  varies  for  each  material.  Above  this  temperature,  all  ferroelectric  crystals  are 
nonpolar,  or  in  a  paraelectric  state.  Ferroelectrics  under  the  Curie  temperature  exhibit  multi¬ 
stable  ground  states.  Therefore,  they  can  change  polarization  directions  through  the  application 
of  an  external  electric  field. 
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Figure  1 .  Schematic  of  the  polarization  phenomenon  in  a  perovskite 
crystal  cell  (3). 

Not  all  crystals  in  a  structure  share  the  same  polarization  directions.  Regions  of  uniform  dipole 
moments  are  referred  to  as  domains,  with  the  areas  between  them  being  domain  walls.  A  wall 
separating  two  anti-parallel  domains  is  a  1 80-degree  wall,  while  a  wall  separating  perpendicularly 
orientated  dipoles  is  a  90-degree  wall  (both  shown  in  figure  2).  The  total  polarization  of  a  ferro¬ 
electric  structure  is  determined  when  the  volume  of  the  domains  is  summed  in  each  direction. 
Ferroelectric  domains  are  formed  to  minimize  the  electrostatic  energy  of  depolarizing  fields  caused 
by  surface  charges  ( 4 ).  The  switching  of  polarization  is  conducted  through  nucleation  and  domain- 
wall  motion.  The  switching  time  is  an  important  property  for  many  applications  and  is  limited  by 
these  two  factors. 


Figure  2.  Polarization  domains  and  domain  walls  in  an 
inhomogeneously  polarized  material  (3). 


Figure  3  shows  that  the  net  polarization  versus  applied  electric  field  produces  a  hysteresis  loop. 

As  an  electric  field  is  applied  to  a  nonpolar  ferroelectric,  domains  begin  to  switch  until  all  domains 
point  in  the  same  direction  (A).  After  the  maximum  polarization  is  reached,  the  electric  field 
decreases  and  some  of  the  domains  switch  back.  At  a  zero  electric  field,  more  domains  still  have  a 
positive  polarization  (B).  The  value  at  this  point  is  the  remnant  polarization.  As  a  negative  field  is 
applied,  more  domains  become  negative  until  the  ferroelectric  is  completely  nonpolar  (C).  The 
field  required  to  reach  a  zero  polarization  is  the  coercive  field.  Finally,  the  field  becomes  more 
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negative  until  a  maximum  negative  polarization  is  reached  (D).  The  process  is  then  reversed  when 
we  switch  back  to  a  positive  electric  field. 


Figure  3.  Hysteresis  loop  common  to  ferroelectric  materials  (3). 

As  a  ferroelectric  switches  back  and  forth  between  the  two  polarized  states,  the  maximum  and 
remnant  polarizations  begin  to  decrease.  This  loss  in  switchable  permanent  polarization  is  ferro¬ 
electric  “fatigue”.  Besides  switching  time  and  ferroelectric  fatigue,  a  third  factor  affecting  the  use 
of  ferroelectrics  is  dielectric  breakdown.  This  occurs  when  an  insulating  material,  such  as  a  ferro¬ 
electric,  is  forced  to  conduct  electricity  in  large  electric  fields.  Dielectric  strength  measures  the 
maximum  electric  field  that  a  material  can  experience  before  dielectric  breakdown  occurs.  The 
strength  is  measured  in  volts  per  unit  length  of  thickness. 


3.  Background  on  Toroidal  Ordering  of  Nanodots 


The  concept  of  toroidal  ordering  of  ferroelectric  nanodots  began  in  1946  with  the  work  of  Charles 
Kittel  with  ferromagnetic  domains  in  small  particles  (5).  His  theory  was  based  on  the  observation 
that  the  demagnetized  state  is  the  stable  state  in  large  ferromagnetic  crystals  in  the  absence  of  an 
applied  magnetic  field.  To  achieve  a  demagnetized  state,  a  complex  ordering  of  the  domains  is 
required.  The  ground  state  can  be  determined  if  the  free  energy  equation  is  optimized  in  terms  of 
the  domain  ordering.  The  free  energy  of  the  structure  is  given  by  the  equation 

F  =  F  +  F  +  F 

1  1  W  ' *  1  M  '  1  a 

The  first  term,  Fw ,  represents  the  surface  energy  between  domains  and  is  calculated  by 

Fw  =<?WS 
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where  aw  is  the  surface  energy  per  unit  area  and  depends  on  the  orientation  of  the  boundary  plane 

and  the  difference  in  spin  direction  across  the  boundary.  The  term,  S,  is  the  total  area  of  the 
domain  boundaries. 

The  second  term  in  the  free  energy  equation,  FM ,  is  the  magnetic  field  energy  of  the  configuration 
or  the  work  required  to  distribute  the  effective  magnetic  poles.  The  energy  is  calculated  by 

where  M  is  the  magnetic  moment  per  unit  volume  (magnetization)  and  H  is  the  resultant  magnetic 
field.  The  integral  is  taken  over  the  volume  of  the  specimen. 

The  final  term  in  the  free  energy  equation,  Fa ,  represents  the  anisotropy  energy  of  spin  orienta¬ 
tion.  This  refers  to  the  energy  associated  with  a  preferred  axis  of  magnetization.  It  is  calculated 
by 

Fa  =  Pa  Va 

where  pa  is  the  anisotropy  energy  density  of  domains  oriented  in  directions  other  than  the  direction 
of  the  easiest  magnetization  and  Va  is  the  volume  of  those  domains. 


Two  cases  for  small  ferromagnetic  nanodots  were  examined  with  the  free  energy  equation.  These 
cases  were  for  a  particle  with  one  domain  and  for  a  particle  with  toroidally  ordered  domains.  The 
first  case  had  only  magnetic  energy  since  the  domain  was  in  the  direction  of  preferred  magnetiza¬ 
tion  and  there  were  no  domain  wall  surfaces.  Therefore,  the  energy  could  be  described  as 


F=  —  M/I1 
3 


The  energy  for  case  two  involved  anisotropy  energy  and  domain  wall  surface  energy  but  had  no 
net  magnetization.  Therefore,  the  energy  could  be  described  as 

F  =  (Twl42L2  +^L3 
2 


Through  this  process,  it  was  discovered  that  the  toroidal  ordering  results  in  the  minimum  energy 
for  L  >  15  nm  using  empirical  estimations  for  the  other  variables. 


4.  Literature  Review 


The  field  remained  relatively  unexplored  over  the  next  30  years  until  Zheludev  et  al.  ( 6)  first 
reported  such  structures.  However,  their  work  could  never  be  reproduced  and  therefore,  the  study  of 
toroidal  ordering  was  again  abandoned.  A  collection  of  papers  between  1984  and  1985  extended 
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Kittle’s  work.  These  papers  include  Ginzburg  et  al.  (7),  Sannikov  (8),  and  Sannikov  and  Zheludev 
(9).  All  these  papers  confirm  the  existence  of  toroidal  moments  in  magnetism.  A  decade  later, 
Gorbatsevich  and  Kopaev  (10)  and  Sannikov  (11)  further  extended  the  work.  However,  the  work 
still  remained  relatively  unnoticed  until  Naumov  et  al.  (2)  performed  ab  initro  simulations  on  freely 
standing  nano-particles  of  lead  zirconate  titanate.  The  simulations  showed  that  phase 
transformations  do  occur  in  zero-dimensional  particles  (nanodots)  and  that  the  phase  transformations 
lead  to  spontaneous  toroidal  moments  under  a  critical  temperature  of  approximately  600K.  These 
toroidal  moments  occur  to  minimize  the  depolarizing  field.  Naumov  et  al.  (2)  found  that  two  stable 
ground  states  exist  for  nanodots  with  a  diameter  larger  than  3.2  nm  and  the  states  can  be  switched  via 
a  magnetic  field  rather  than  a  static  electric  field  as  used  for  bulk  ferroelectrics. 

The  work  was  continued  by  papers  such  as  Ponomareva  et  al.  (11,  12)  and  Prosandeev  et  al.  (13). 
Ponomareva  et  al.  published  two  papers  in  2005  proposing  models  to  examine  the  depolarizing 
fields  (4)  and  effects  of  electrical  and  mechanical  boundary  conditions  on  ferroelectric  nanoparti¬ 
cles  (12).  The  atomistic  models  examine  the  depolarizing  field  and  energy  for  ferroelectric  nano¬ 
structures.  They  found  the  depolarizing  field  and  energy  by  calculating  the  difference  between 
dipole-dipole  interactions  for  ferroelectric  nanostructures  during  perfect  open-circuit  and  short- 
circuit  electrical  boundary  conditions.  An  open-circuit  condition  exhibits  a  maximum  depolarizing 
field,  while  a  short-circuit  condition  completely  screens  the  effect  of  the  depolarizing  field.  The 
energies  calculated  via  their  model  were  independent  of  dot  size  and  agreed  exactly  with  the 
established  continuum  approach.  Through  their  simulations,  Ponomareva  et  al.  found  that  toroidal 
moments  occur  when  the  screening  of  the  depolarizing  field  is  less  than  88%  for  a  freely  standing 
lead  zirconate  titanate  nanostructure.  After  applying  tensile  and  compressive  strains  of  2.65%, 
they  discovered  that  toroidal  moments  occur  under  tensile  strain  with  a  screening  of  less  than  88% 
but  not  under  compressive  strain.  These  results  are  illustrated  in  figure  4. 

Prosandeev  et  al.  (13)  proposed  a  switching  mechanism  that  uses  inhomogeneous  static  electric 
fields  to  alter  the  toroidal  moment  of  the  nanodots  without  affecting  the  surrounding  nanodots. 

They  did  this  by  placing  two  opposite  dipolar  sources  at  each  end  of  the  dot.  This  setup  prevents 
any  net  polarization  that  could  affect  nearby  nanodots.  By  rotating  the  source,  they  could  switch 
the  toroidal  moment.  The  ability  to  switch  the  direction  of  the  toroidal  moment  is  vital  for  the 
integration  of  ferroelectric  nanodots  (FENs)  into  memory  devices  and  other  applications. 
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Figure  4.  (a)  Freely  standing  nanodot  (b)  FEN  under  tensile  strain  (c)  FEN  under  compressive  strain  (72). 


5.  Advantages  of  Ferroelectric  Nanodots 


An  obstacle  for  implementing  FENs  into  NVFRAM  is  the  understanding  of  the  phase  transforma¬ 
tions  on  the  nanoscale.  In  certain  conditions,  the  polarization  of  nanodots  is  oriented  toroidally 
rather  than  rectilinearly.  The  local  dipoles  are  directed  in  a  suitable  direction  to  equilibrate  the 
depolarizing  field  (2).  The  resultant  orientation  produces  a  toroidal  moment.  Various  properties  of 
the  toroidally  ordered  nanodots  are  conducive  for  use  in  NVFRAM.  The  phase  is  bistable,  meaning 
that  two  minimum  energy  states  exist.  These  states  can  be  easily  translated  to  “0”  or  “1”  in  Boolean 
logic.  Ground  states  can  be  switched  with  the  use  of  a  magnetic  field  that  does  not  require  contact 
with  electrodes,  which  is  challenging  to  create  on  the  nanoscale  (2).  The  local  dipole  orientation 
does  not  result  in  a  net  polarization;  therefore,  neighboring  nanodots  will  not  be  affected  when  one 
nanodot  switches  polarization.  All  these  properties  make  toroidally  ordered  FENs  appealing  to 
NVFRAM. 
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6.  Future  Work 


The  future  of  toroidally  ordered  ferroelectric  nanodots  in  memory  devices  is  bright.  In  principle, 
combining  the  properties  shared  by  all  ferroelectrics  with  those  unique  to  toroidal  ordering  can  lead 
to  smaller,  faster,  and  more  reliable  memory  materials.  However,  more  research  must  be  performed 
on  all  aspects  of  their  properties  and  behaviors  before  ferroelectric  nanodots  can  be  implemented. 
Among  many  mechanisms  that  remain  unstudied,  the  interacting  relationships  among  atomic 
structure,  electromagnetic  fields,  and  mechanical  strains  are  not  known.  For  example,  the  role  of 
mechanical  strains  and  thermal  expansion  on  the  stability  of  the  toroids  is  unknown.  Most 
importantly,  it  is  still  unclear  if  toroidal  ordering  has  ever  occurred  experimentally  since  concepts 
have  been  posited  mainly  through  theoretical  models  and  calculations.  Therefore,  experimentation 
must  accompany  the  theoretical  calculations  before  ferroelectric  nanodots  are  seriously  considered 
for  NVFRAM. 
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Appendix  A.  Derivation  of  Ewald  Sum  Equation 


The  following  mathematical  derivation  was  an  assignment  during  the  internship.  The  details  of  the 
method  were  uncovered  to  help  us  understand  basic  pedagogy  of  modeling  ferroelectric  systems 
from  the  atomistic  level. 


The  Ewald  sum  is  a  means  to  calculate  the  electrostatic  interactions  of  periodic  systems.  Its  contri¬ 
bution  to  the  total  energy  enables  the  formation  of  long  range  order  and  is  the  physical  underpinning 
of  ferroelectric  polarization.  The  method  involves  superimposing  a  spherical  Gaussian  cloud  on 
each  ion  in  real  space.  These  spherical  clouds  have  a  charge  opposite  the  ion.  A  spherical  Gaussian 
cloud  is  also  superimposed  on  each  ion  in  reciprocal  space.  These  clouds  have  the  same  charge  as 
the  ion  they  surround.  The  Ewald  sum  calculation  involves  using  Poisson’s  equation  and  Fourier 
transforms.  The  Ewald  sum  equation  is  comprised  of  a  term  to  calculate  short-range  interactions  in 
real  space,  a  term  to  calculate  long-range  interactions  in  reciprocal  space,  and  a  correction  term  for 
self  energy. 

General  Equation  for  Electrostatic  Energy: 


t'-lzzs 

^  i  j  n 


ri  -  ( rj  +  nL) 


where  1/2  accounts  for  double  counting  of  atom  interactions 

qi  and  qj  are  the  charges  on  atoms  i  and  j 
ri  and  rj  are  the  coordinates  of  atoms  i  and  j  in  the  unit  cell 
nL  is  the  displacement  vector  for  the  different  image  cells 


where 


<f)(n)  =  Electric  Potential  at  ri 
V, 


/ 

n  rf  -  (r.  +  nL) 


Charge  Density: 

p(r)  =  Point  Charge  -  Gaussian  Charge  +  Gaussian  Charge 
Density  of  Point  Charge:  p(r)  =  Y^q.Sir-r) 

i 

where  S(r  -  rt )  is  the  Dirac  delta  function 
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Density  of  Gaussian  Charge:  p{r )  =  VV  — 

nj 


\  3/2 


- a(r-rf ) 


where  a  is  a  parameter  affecting  the  Gaussian  width  and  affects  computational  efficiency 

P(r)  =  YJ< 


( 

/  „  \ 

3/2  A 

/  \  3/2 

<cti 

1 

1 

1 

£ 

>i 

1 

+  qt 

— -  e~a(r~ri)2 

V 

) 

Kn) 

J 


V  Y~ 

Short  Range  Contributions  Long  Range  Contributions 

Electric  Potential  (Sh 

<j> short  (r)  =  Potential  from  Point  Charge  -  Potential  from  Gaussian  Charge 


Potential  from  Point  Charge:  <f>short  P  (r)  =  ^ 


Potential  from  Gaussian  Charge: 
Poisson’s  Equation:  V2^(r)  =  4 np{r) 


<h 


(r-r,) 


p(r)  =  qt 


/  \3/2 

'  a  ' 


j 


- a{r-rt ) 


1  d2r^(r) 
r  dr 2 

d  r</){r ) 
dr 


=  4  np{r) 


f  4 nr  p(r)  dr 

Joo 

/•OO 

-Jr  2 nqt 


r  \3/2 

“I  e^dr2 
W 


/  \  3/2 

a 


■2nq\  —  r  e-arldr2 

U )  Jr 


■2  qt 


/  \  1/2 
'  a  N 


/  \l/2 

'  a  ' 


r(f>(r)  =  2qt  [e^  dr 

W  Jo 

(j){r)  =  —erf  (-Ja  r) 
r 

where  erf  is  the  error  function  and  r=(r-ri) 

2L»  =  ZI| 


Vi 


■(ri+nL)  |  |  r  —  (ri+  nL)  | 


erf(ja(r-(r  +nL ))) 
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<h 


-  erfc{\[a{r  -  (ji  +  nL))) 


l'n\r~(ri+nL)\ 

where  erfc  is  the  complimentary  error  function 


Electrostatic  Energy  (Short  Range): 


1 


‘111 


J 


Z'ty  n  \r  -  (r .  +  nL) 


Correction  Term  for  Self  Interaction: 


erfc[-Ja  (rt  -  (r^.  +  nL ))) 


<f>{r)  =  —  erf  (-Ja  r ) 


0(r  =  0)  =  qt 


f  2  ' 


vV^  J 


\e'ar  yja 


By  taking  the  limit  as  r  approaches  zero  (using  L’Hopital), 

</>self(ri)  =  2<li 

U self  =  i^elf  (r,  ) 


/  X 1/2 

'  a 
\n) 


1 


Uself  = 


a  2 

\K  P 


2>/ 


Electric  Potential  (Reciprocal  Space): 

Poisson’s  Eq:  -  V2^(r)  =  4 np(r) 
Poisson’s  Eq  after  Fourier  Transform: 

</>(k)  =  jrP(k) 

K 


k2(j){k )  =  \7zp{k) 


kylr 


P(k)  = 


i  lp(r)e*dr 

J 


V 


e~ar  e  -,k'  dr 


if 

V 


J  n 
3 

aV 

- 

y  all  space 


/  2 

a(r  H - ) 


Jr 
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Side  note: 


f 

7 

ikr') 

r  an 

r  H - 

= 

r  H - 

V 

a  J 

l  2a ; 

+  ¥ 


4  a" 


V 


(-Y 

Kn) 


f 


ik  X  kz 

r+ -  + - - 

2  a)  4  a2 


Xr 


q_L 

v 


r~'i  =* 


Since  the  bounds  of  the  integral  are  -  oo,  oo , 

<•00  2 

f  e  ** 

J-GO 

n 

{erf  (oo)  -  erf  (— oo)) 


(  ik 
r  H - 

V  2  aj 


can  be  replaced  by  r2 


a 

\7Z  j 


,  4a 


f_L 

V 


/  \  —  -k 

a  p  -7- 
—  e 4 


2  fa 


f 

V 


^  af 


- k 
i  4  a 


\7T  J 


1  — 

px)=jYj^je'a 


</>(k)  =  fTjHqje 


-k 
4  a 


(j)(r)  =  YJ<t>(k)e 

k 

4k 


ikr 


V  k  j  & 


=  2ffyyy 

lons  v  r  tj  x 


ik*(ri-rj) 


-k 
7  4  a 


In  ref  12  —  is  omitted 

K 

because  of  choice  of  units. 
This  term  does  not  appear 
in  any  other  derivation 
either  and  is  attributed  to 
the  choice  between  SI  and 
CGS. 


Ewald  Sum: 


u  =  ushort+u. 


long 


■U 


self 


in 


m 


i  ij 


erfc(fa(r  -  (rj  +  nL)))+  ^  X Z  Z 


r 


■  (jj  +  nL) 


M 


V 


•~JJ  e  4 a  eik<r< -0) 
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